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Abstract 

A study on the effects of optical gain nonuniformly distributed in one-dimensional random sys- 
tems is presented. It is demonstrated numerically that even without gain saturation and mode 
competition, the spatial nonuniformity of gain can cause dramatic and complicated changes to 
lasing modes. Lasing modes are decomposed in terms of the quasi modes of the passive system 
to monitor the changes. As the gain distribution changes gradually from uniform to nonuniform, 
the amount of mode mixing increases. Furthermore, we investigate new lasing modes created 
by nonuniform gain distributions. We find that new lasing modes may disappear together with 
existing lasing modes, thereby causing fluctuations in the local density of lasing states. 
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I. INTRODUCTION 



Lasing modes in random media behave quite differently depending on the scattering 
characteristics of the media PQ. In the strongly scattering regime, lasing modes have a 
nearly one-to-one correspondence with the localized modes of the passive system j2j [3] . Due 
to small mode volume, different localized modes may be selected for lasing through local 
pumping of the random medium (2], H]. The nature of lasing modes in weakly scattering 
open random systems (e.g., [5j EJ [7]) is still under discussion [8]. In systems which are 
diffusive on average, prelocalized modes may serve as lasing modes In general, however, 
the quasi modes of weakly scattering systems are very leaky. Hence, they exhibit a large 
amount of spatial and spectral overlap. For inhomogeneous dielectric systems with uniform 
gain distributions, even linear contributions from gain induced polarization bring about 
a coupling between quasi modes of the passive system [10J. Thus, lasing modes may be 
modified versions of the corresponding quasi modes. However, Vanneste et al. found that 
when considering uniformly distributed gain, the first lasing mode appears to correspond to 
a single quasi mode [11J. The study was done near the threshold pumping rate and nonlinear 
effects did not modify the modes significantly. Far above threshold, it was found that lasing 
modes consist of a collection of constant flux states [12]. Mode mixing in this regime is 
largely determined by nonlinear effects from gain saturation. 

Remaining near threshold, pumping a local spatial region, and including absorption out- 
side the pumped region found lasing modes to differ significantly from the quasi modes of 
the passive system [13]. This change is attributed to a reduction of the effective system 
size. More surprisingly, recent experiments [131 IT5] and numerical studies [16] showed the 
spatial characteristics of lasing modes change significantly by local pumping even without 
absorption in the unpumped region. It is unclear how the lasing modes are changed in 
this case by local pumping. In this paper, we carry out a detailed study of random lasing 
modes in a weakly scattering system with a nonuniform spatial distribution of linear gain. 
Mode competition depends strongly on the gain material properties, e.g., homogeneous vs. 
inhomogeneous broadening of the gain spectrum. Ignoring gain saturation (usually respon- 
sible for mode competition) and absorption, we find that spatial nonuniformity of linear 
gain alone can cause mode mixing. We decompose lasing modes in terms of quasi modes 
and find them to be a superposition of quasi modes close in frequency. The more the gain 
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distribution deviates from being uniform, more quasi modes contribute to a lasing mode. 

Furthermore, still considering linear gain and no absorption outside the gain region, we 
find that some modes stop lasing no matter how high the gain is. We investigate how the 
lasing modes disappear and further investigate the properties of new lasing modes (TTj that 
appear. The new lasing modes typically exist for specific distributions of gain and disappear 
as the distribution is further altered. They appear at various frequencies for several different 
gain distributions. 

In Section [XXJ, we describe the numerical methods used to study the lasing modes of a 
one-dimensional (ID) random dielectric structure. The model of gain and a scheme for 
decomposing the lasing modes in terms of quasi modes is presented. A method to separate 
traveling wave and standing wave components from the total electric field is introduced. 



The results of our simulations are presented and discussed in Section III Our conclusions 



concerning the effects of nonuniform gain on lasing modes are given in Section IV 



II. NUMERICAL METHOD 

The ID random system considered here is composed of N layers. Dielectric layers with 
index of refraction n > 1 alternate with air gaps (n = 1) resulting in a spatially modulated 
index of refraction n(x). The system is randomized by specifying different thicknesses for 
each of the layers as = (^1,2) (1 + vO where (di) and (d 2 ) are the average thicknesses 
of the layers, < r\ < 1 represents the degree of randomness, and ( is a random number in 

(-1,1)- 

A numerical method based on the transfer matrix is used to calculate both the quasi 
modes and the lasing modes in a random structure. Electric fields on the left (right) side of 
the structure p (q^) and q (pn) travel toward and away from the structure, respectively. 
Propagation through the structure is calculated via the 2x2 matrix M 



(1) 



The boundary conditions for outgoing fields only are p = q^ = 0, requiring M 22 = 0. 

For structures without gain, wavevectors must be complex in order to satisfy the 
boundary conditions. The field inside the structure is represented by p(x) exp[in(x)kx] + 
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q(x) exp[—in(x)kx], where k = k + ik{ is the complex wavevector and x is the spatial coor- 
dinate. For outgoing-only boundary conditions (M 22 = 0), ki < 0. The resulting field dis- 
tributions associated with the solutions for these boundary conditions are the quasi modes 
of the passive system. 

Linear gain is simulated by appending an imaginary part to the dielectric function e(x) = 
e r (x) + iei(x), where e r (x) = n 2 (x). This approximation is only valid at or below threshold 
[18] . In Appendix A the complex index of refraction is calculated as n(x) = y~e{ 
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n r (x) + irii, where Hi < 0. We consider to be constant everywhere within the random 
system. This yields a gain length l g = 1/ki = 1/riik (k = 2tt/X is the vacuum wavenumber 
of a lasing mode) which is the same in the dielectric layers and the air gaps. Note that 
this gain length is not to be confused with the length of the spatial gain region, which is 
described below. The real part of the index of refraction is modified by the imaginary part 
) = yjn 2 (x) + rz|. A real wavevector k = 2ir/\ describes the lasing frequency. The 
field inside the structure is now represented by p(x) exp[in(x)kx] + q(x) exp[— in{x)kx\. The 
frequencies and thresholds are located by determining which values of k and n*, respectively, 
satisfy M 22 = 0. 

Nonuniform gain is introduced through an envelope function /^(x) multiplying rij. The 
envelope considered here is the step function Je{x) = H(—x + Iq), where x = is the left 
edge of the random structure and x = Iq is the location of the right edge of the gain region. 
Ig may be chosen as any value between and L. 

The solutions of the system are given by the points at which the complex transfer matrix 
element M 22 = 0. Where Re[M 22 ] = or Im[M 22 ] = 0, "zero lines" are formed in the plane of 
either complex k (passive case) or (k, rii) (active case). The crossing of a real and imaginary 
zero line results in M 22 = at that location, thus revealing a solution. We visualize these 
zero lines by plotting log 10 |ReM 22 | and log 10 |ImM 22 | to enhance the regions near M 22 = 
and using various image processing techniques to enhance the contrast. 

The solutions are pinpointed more precisely by using the Secant method. Locations of 
minima of |M 22 | 2 and a random value located closely to these minima locations are used as 
the first two inputs to the Secant method. Once a solution converges or |M 22 | < 1CT 12 , a 
solution is considered found. This method has proved extremely adept at finding genuine 
solutions when suitable initial guesses are provided. 

Verification of these solutions is provided by the phase of M 22 , calculated as 9 = 
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atan2(ImM22, Reil^)- Locations of vanishing M22 give rise to phase singularities since 
both the real and imaginary parts of M 2 2 vanish. The phase change around a path sur- 
rounding a singularity in units of 2ir is referred to as topological charge [T51 |2"U] . In the 
cases studied here, the charge is +1 for phase increasing in the clockwise direction and — 1 
for phase increasing in the counterclockwise direction. 

Once a solution is found, the complex spatial field distribution may be calculated. Quasi 
modes ip(x) are calculated in the passive case and lasing modes ty(x) are calculated in the 
active case. In order to determine whether or not mode mixing occurs (and if so, to what 
degree) in the case of nonuniform gain, the lasing modes are decomposed in terms of the 
quasi modes of the passive system. It was found [2H [22] that any spatial function defined 
inside an open system of length L [we consider the lasing modes ^>(x) here] can be expressed 
as 

^(X) = J2 a m^Pm(x), (2) 

m 

where ip m ( set of quasi modes characterized by the complex wavevectors k m . The 

coefficients calculated by 

a m =i / ^(x)ip m (x) + $(x)ip m (x) dx 
Jo 1 J 

+ z[*(0)^ m (0) + *(L)^ m ( J L)], (3) 

where 

i> m {x) = -ik m n 2 {x)i) m (x) (4a) 
= -ik [n r (x) + in i f E (x)] 2 ty(x). (4b) 

The normalization condition is 



1 =i / 2^{x)^j{x)dx + i [ip 2 {0) + ip 2 {L)] 
Jo 

=i I 2^(x)^(x)dx + i [^ 2 (0) + ^ 2 (L)] . (5) 
Jo 

An advantage of this decomposition method is that a calculation over an infinite system 
has been reduced to a calculation over a finite system. Error checking is done by using 
the coefficients found in Eq. ^ to reconstruct the lasing mode intensity distribution with 
Eq. ([2]) yielding R(x) = ^ m a m il) m (x). We define a reconstruction error E R to monitor the 
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accuracy of the decomposition: 



_ f\*(x)-R(x)\ 2 dx 



In general, the field distribution of a mode in a leaky system consists of a traveling 
wave component and a standing wave component. In Appendix [B] we introduce a method 
to separate the two. This method applies to quasi modes as well as lasing modes; we shall 
consider lasing modes here. At every spatial location x, the right-going complex field ty( R \x) 
and left-going complex field ^^(x) are compared. The field with the smaller amplitude is 
used for the standing wave and the remainder for the traveling wave. If |^^(a;)| < |^^(x)| 
then the standing wave component ^^ s \x) and traveling wave component ^^ T \x) are 

q?( s ) ( x ) = q?( R ) ( x ) + (x)}* (7a) 

= q?M( x ) _ (7b) 



Further physical insight on lasing mode formation and disappearance, as well as new 
lasing mode appearance is provided by a mapping of an "effective potential" dictated by the 
random structure. Local regions of the random medium reflect light at certain frequencies 
but are transparent to others [23J. The response of a structure to a field with frequency 
uj = ck can be calculated via a wavelet transformation of the real part of the dielectric 
function e r {x) = n 2 (x) |24J. The relationship between the local spatial frequency q res and 
the optical wavevector k is approximately q res = 2k in weakly scattering structures. The 
Morlet wavelet \ is expressed as 

x ( x ' ~ x \ = 7r ~ 1 4 c W«'-»)/« c -fa'-») 3 /2« 2 ; ( 8 ) 



with nondimensional frequency uoq and Gaussian envelope width s [25]. With c^o fixed, 
stretching the wavelet through s changes the effective frequency. Wavelets with varying 
widths are translated along the spatial axis to obtain the transformation 

W(x,q res ) = J e r {x') X * (^^) dx'. (9) 

where 



ft- = + (10) 
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FIG. 1: (Color online) Effective potential (wavelet power spectrum) |H^(x)| 2 of the dielectric 
function e r {x) as a function of position x and wavelength A. Regions of high power indicate potential 
barriers and regions of low power indicate potential wells where intensities are typically trapped. 
The black lines on the top represent the spatial distribution of dielectric constant e r (x) = n 2 (x). 

The wavelet power spectrum \W(x,q res )\ 2 is interpreted as an effective potential. Regions 
of high power indicate potential barriers and regions of low power indicate potential wells 
for light frequency u> = q res c/2. 

III. RESULTS 

A random system of N = 161 layers is examined in the following as an example of a 
random ID weakly scattering system. The indices of refraction of the dielectric layers are 
Hi = 1-05 and the air gaps ri2 = 1. The average thicknesses are (di) = 100 nm and (o^) = 200 
nm giving a total average length of (L) = 24100 nm. The grid origin is set at x = and the 
length of the random structure L is normalized to (L). The degree of randomness is set to 
rj = 0.9 and the index of refraction outside the random media is n = 1. The localization 
length £ is calculated from the dependence of transmission T of an ensemble of random 
systems of different lengths as = — d(lnT) jdL. The above parameters ensure that the 
localization length is nearly constant at 200 fim < £ < 240 /xm over the wavelength range 
500 nm < A < 750 nm. With £ 3> L, the system is in the ballistic regime. 
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FIG. 2: (Color online) A mapping of the phase 9 of M22 for the passive ID random system without 
gain. The topological charge of all quasi modes seen here is —1. Modes are enumerated from left 
to right. Quasi mode 27 is encircled in black. 

Figure [T] shows the effective potential of the structure within the wavelength range of 
interest via a wavelet transformation. We use a nondimensional frequency of ujq = 6 [26J and 
a spatial sampling step of Ax = 2 nm. The power spectrum |PF(a;)| 2 reveals the landscape 
of the effective potential dictated by the locations and thicknesses of the dielectric layers. 

Figure [2] is a phase map of M 2 2 in the passive case (without gain). The phase singularities 
mark the quasi modes' k values and are indicated by phase changes from — 71 to tt along any 
lines passing through. The topological charge of all quasi modes is —1. Adjacent modes 
are formed by real and imaginary zero lines of M22 that are not connected to one another. 
We calculated M 2 2 for increasingly large \k{\ values until machine precision was reached and 
no additional modes appeared. As previously found P~BJ, mode frequency spacing is fairly 
constant in the ballistic regime. The nearly equal spacing of phase singularities in Fig. [2] 
attests to this. 

Most quasi modes have similar decay rates except for a few which have much larger decay 
rates. Modes are enumerated here starting with the lowest frequency mode in our wavelength 
range of interest. Mode 1 has a wavelength of 748 nm and mode 33 has a wavelength of 
502 nm. Most quasi modes have k{ values around —0.1 [imr 1 . But a few have much larger 
decay rates, such as mode 27 at A = 532 nm which has fcj = —1.03 /im -1 (encircled in black 
in Fig. [2]). Figure [3] shows the intensity of mode 27 to be concentrated on one side of the 
open structure. We observe that it bears similarity to "doorway states" common to open 



8 



I 

1 

1 

y. 



5000 10000 15000 20000 

* (nm) 

FIG. 3: Normalized intensity ^(rc)! 2 of a leaky mode-quasi mode 27 (k = 11.8 — il.03 ^m^ 1 ) of 
the passive random ID structure. The intensity is peaked at the left boundary of the structure, 
similar to doorway states |27j . 

quantum systems [27] . Doorway states are concentrated around the boundary of a system 
and strongly couple to the continuum of states outside the structure. Therefore, they have 
much larger decay rates. 

For the case of uniform gain, only the lasing modes with large thresholds change sig- 
nificantly from the quasi modes of the passive system. Finding the corresponding quasi 
modes for lasing modes with large thresholds is challenging due to changes caused by the 
addition of a large amounts of gain. Thus, we neglect them in the following comparisons. 
However, there is a clear one-to-one correspondence with quasi modes for the remaining 
lasing modes. The average percent difference between quasi mode frequencies and lasing 
mode frequencies is 0.026%. The average percent difference between quasi mode decay 
rates ki and lasing thresholds (multiplied by k for comparison) riik is 4.3%. The normal- 
ized intensities of the quasi modes Iq(x) = \tp(x)\ 2 and lasing modes = \^(%)\ 2 are 
also compared. The spatially averaged percent difference between each pair of modes is 
calculated as (2/L) J {\Iq(x) — Il(x)\/[Iq(x) + Il(x)]) dx x 100. The averaged difference 
between intensities of the 3 lasing modes with the largest thresholds (of the large threshold 
modes not neglected) compared to the quasi modes is 68% while the remaining pairs average 
a difference of 4.0%. 
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FIG. 4: (Color online) Frequencies and thresholds (fc, n,) of the lasing modes of the ID random 
structure with gain. Lasing modes 1, 17, 18, and 33 are explicitly marked. The gain region length 
Iq reduces from uniform gain [Iq = L) to nonuniform gain Iq < L. The color indicates the value 
of Iq (units of nm) decremented along the layer interfaces. Due to the random thicknesses of the 
layers, the Iq decrements are unequal. Hence, the reason for the unequal spacing of the color code. 

A. Nonuniform Gain Effects on Lasing Mode Frequency, Threshold, and Intensity 
Distribution 

Figure [1] maps the (k, rii) values of lasing modes as nonuniform gain is introduced by 
reducing the gain region length from Iq = L. In this weakly scattering system the inten- 
sity distributions of modes are spatially overlapping. This results in a repulsion of mode 
frequencies [28J . As the size of the gain region changes, the envelopes of the intensity distri- 
butions change, but for most modes rii is small enough to leave the optical index landscape 
unchanged. Thus, the modes continue to spatially overlap as the size of the gain region 
changes and their frequencies remain roughly the same as in the uniform gain case. Similar 
behavior of lasing mode frequencies can be seen as the gain region length is varied in a 
simpler cavity with uniform index. Thus, the robustness of frequency is not due to inhomo- 
geneity in the spatial dielectric function. However, the threshold values of the lasing modes 
change as Iq decreases. Due to the limited spatial region of amplification, the thresholds 
increase. The increase of rii due to the change of threshold, though considerable, is not large 
enough to significantly impact the lasing frequencies as evidenced by the small change of 
frequencies as Iq decreases. 
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FIG. 5: (Color online) Normalized intensity of lasing mode 17 with uniform gain Iq = 24100 nm 
(red solid lines) and nonuniform gain Iq = 14284 nm (black dashed lines). Gain is only located 
in the region < x < Iq. (a) Total intensity |^(x)| 2 , (b) traveling wave intensity ^^(a:)! 2 , and 
(c) standing wave intensity |^^(a;)| 2 . Nonuniform gain significantly changes the spatial intensity 
envelope as well as the standing wave and traveling wave components. 

The intensity distributions of the lasing modes also change considerably as Iq is reduced. 
Normalized spatial intensity distributions are given by |\l/(a;)| 2 after ty(x) has been normal- 
ized according to Eq. (|5|. The intensities are sampled with a spatial step of Ax = 2 nm. 
With uniform gain (Iq = 24100 nm), the intensity of lasing mode 17 (A = 598 nm) in Fig. 
[5](a) increases toward the gain boundaries due to weak scattering and strong amplification. 
When the gain boundary is changed to Iq = 14284 nm, the envelope of the spatial intensity 
distribution changes dramatically. The intensity increases more rapidly toward the bound- 
aries of the gain region and stays nearly constant outside the gain region but still inside the 
structure. This change can be understood as |rij| inside the gain region causes the intensity 
to become larger, while outside the gain region = and the wavevector is real. 

To monitor the change in the trapped component of the intensity, ^>(x) is separated into 
a traveling wave and a standing wave component via Eq. ^ (see Appendix [Bj. Figures 
[5^b) and (c) show the traveling wave and standing wave components of lasing mode 17, 
respectively. For Iq = L, the intensity increase toward the structure boundaries is caused 
mostly by the growth of the traveling wave. The standing wave part is strongest near 
the center of the system. For Iq = 14284 nm, the standing wave exhibits two peaks, 
one concentrated near the center of the gain region and the other outside the gain region. 



11 




I" 



10000 



)0 20000 



1 1 



" 



5000 



15 



FIG. 6: (Color online) Standing/traveling wave ratio Ast(x) (red solid lines) of lasing mode 
17 with uniform gain Iq = L (a) and nonuniform gain Iq = 14284 nm (b). Gain is located in 
the regions left of the vertical black solid line. We term the location at which the ratio diverges 
(Ast{x) — ► oo), the standing wave center (SWC). The potential profile |VF(x)| 2 (black dashed 
lines) of the dielectric function at the wavelength of mode 17 is overlaid in both (a) and (b). 

However, the standing wave intensity outside the gain region should not be directly compared 
to the standing wave intensity inside the gain region. The total intensity inside the gain 
region increases toward the gain boundary in this weakly scattering system. Thus, the 
amplitude of the field outside the gain region, where there is no amplification, is determined 
by the total field amplitude at the gain boundary. The randomness of the dielectric function 
outside the gain region traps part of the wave which results in a relatively large standing 
wave intensity compared to inside the gain region. However, outside the gain region, there 
is a net flux toward the right boundary of the system meaning the traveling wave intensity 
in this region is large as well. 

The relative strength of the standing wave is calculated through the ratio of standing 
wave amplitude to traveling wave amplitude. The amplitudes are calculated in Appendix |B} 
Depending on whether the prevailing wave is right-going or left-going, the standing/ traveling 
wave ratio is given by 



Results from considering uniform and nonuniform gain for lasing mode 17 are shown in Fig. 




(11) 
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[6j Where the standing wave is largest inside the gain region, |\I/^(a;)| = |0(x) — Il(x)| = 
and the ratio Ast(x) is infinite. The location where the ratio is diverging is the position of 
pure standing wave. Fields are emitted in both directions from this position. The prevailing 
wave to the right of this standing wave center (SWC) is right-going. The prevailing wave to 
the left of this SWC is left-going. The SWC of the lasing mode is located near the center 
of the total system when considering uniform gain in Fig. |6^a). With the size of the gain 
region reduced in Fig. pTb), we see that the SWC of the lasing mode (where Ast(x) — > oo) 
moves to stay within the gain region. Furthermore, note that this mode has a relatively 
small threshold (see Fig. [4]). We have found that in general, modes with low thresholds 
have a SWC near the center of the gain region while high threshold modes have a SWC near 
the edge of the gain region. 

The cause for the small peak of A$t{x) outside the gain region can be found in the 
potential profile of Fig. [TJ A slice of the potential profile ^(x)) 2 at the wavelength of mode 
17 (A = 598 nm) is overlaid on the intensities in Fig. [6] This suggests the standing wave is 
weakly trapped in a potential well around x = 20500 nm [marked by an arrow in Fig. p[b)]. 




B. Mode Mixing 

Lasing modes can be expressed as a superposition of quasi modes of the passive system 
via Eq. ^ for any distribution of gain. Coefficients obtained from the decomposition of 
the lasing modes in terms of the quasi modes by Eq. (|3| offer a clear and quantitative 
way to monitor changes of lasing modes by nonuniform gain. Using Simpson's rule for the 
numerical integrations and a basis consisting of at least 15 quasi modes at both higher and 
lower frequencies than the lasing mode being decomposed, we consistently find Er w 10~ 4 . 

Figure |7^a) shows the decomposition of lasing mode 17 with uniform and nonuniform 
gain. Beginning with the case of uniform gain (Iq = L), the largest contribution to lasing 
mode 17 is from corresponding quasi mode 17. There is a nonzero contribution from other 
quasi modes on the order 10~ 3 . This reflects slight differences between the lasing mode 
profile in the presence of uniform gain and the quasi mode profile [TQl EE] . With the gain 
region length reduced to Iq = 14284 nm, the coefficients \a m \ 2 from quasi modes closer in 
frequency to the lasing modes increase significantly; i.e., quasi modes closer in frequency 
are mixed in. The exceptions are the very leaky quasi modes 7, 14, and 23. Unlike leaky 
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FIG. 7: (Color online) Decomposition of lasing mode 17 in terms of passive quasi modes, (a) 
Decomposition with uniform gain (red crosses) and nonuniform gain (black circles). Leaky quasi 
modes, i.e., modes with large \ki\ such as modes 7, 14, and 23, contribute to lasing modes markedly 
different than the others, (b) Five largest coefficients from the decomposition of lasing mode 17. 
As Iq is reduced, the amount of mode mixing increases dramatically. The reconstruction error Er 
for lasing mode 17 is close to 10 -4 until Iq = 11000 nm then rises to 10 -2 at Iq = 3200 nm. Some 
coefficients are greater than one, which is possible in open systems. 

quasi mode 27 shown in Fig. [3j quasi modes 7, 14, and 23 have intensities which are peaked 
at the right boundary of the structure. It has been observed that when Iq reduces and the 
intensity distribution of lasing mode 17 shifts to the left boundary of the structure, there is 
less overlap with these leaky quasi modes. Thus, the magnitude of the coefficients associated 
with the leaky modes reduces as shown in Fig. (7^a). 

Figure [7]^b) reveals the five largest coefficients \a m \ 2 for lasing mode 17 as Iq is incre- 
mentally reduced along the dielectric interfaces. While the lasing mode remains dominantly 
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composed of its corresponding quasi mode, neighboring quasi modes mix into the lasing 
mode significantly. It was shown [TU] that linear contributions from gain induced polariza- 
tion bring about a coupling between quasi modes of the passive system. This coupling arises 
solely due to the inhomogeneity of the dielectric function, not the openness of the system. 
While this interaction may play a role in mode mixing with uniform gain, the effect is small 
compared to the mode mixing caused by the nonuniformity of the gain. This is clearly 
demonstrated in Fig. [7](b), where the coefficients from quasi modes close in frequency are 
orders of magnitude larger for small Xq than for Iq = L. 

C. Lasing Mode Disappearance and Appearance 

As the size of the gain region reduces we observe that some lasing modes disappear 
and new lasing modes appear. The existence of new lasing modes in the presence of gain 
saturation has been confirmed |17j . This phenomenon is not limited to random media, but 
its occurrence has been observed in a simple ID cavity with a uniform index of refraction. 
New lasing modes, to the best of our knowledge, are always created with larger thresholds 
than the existing lasing modes adjacent in frequency. The disappearance of lasing modes 
is not caused by mode competition for gain because gain saturation is not included in our 
model of linear gain. Disappearance/appearance events occur more frequently for smaller 
values of Ig- New lasing modes appear at frequencies in between the lasing mode frequencies 
of the system with uniform gain. These new modes exist only within small ranges of Ig- We 
also find that the disappearance events exhibit behavioral symmetry (as explained below) 
around particular values of Ig- This disappearance and subsequent reappearance causes a 
fluctuation of the local density of lasing states as Ig changes. 

We examine the progression of one representative event in detail. The gaps in the decom- 
position coefficients for lasing mode 17 in Fig. ^b), in the range 10500 nm < Ig < 14500 
nm, indicate lasing mode 17 does not exist for those distributions of gain. Figure [8] shows the 
real and imaginary zero lines of M22 and their accompanying phase maps for Ig = 14961 nm, 
14553 nm, 14523 nm, 14472 nm, 14284 nm, and 14042 nm. As Iq decreases, the zero lines 
of lasing modes 17 and 18 join as seen in the transition from Fig. |8^a) to (c). This creates 
a new mode solution (marked by a white circle) with a frequency between lasing modes 17 
and 18 and a larger threshold. The existence of a new lasing mode is confirmed by the phase 
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FIG. 8: (Color online) (left) Real (green) and Imaginary (red) zero lines of Mi2- Their crossings 
indicate (k, rii) values of lasing modes, (right) Phase maps of M22. All data are in the ranges (10.3 
fiia < k < 10.8 iivaT 1 ) and (0 > > —0.074) covering lasing modes 17 and 18 for Iq = 14961 
run (a-b), 14553 nm (c-d), 14523 nm (f-g), 14472 nm (h-i), 14284 nm (j-k), and 14042 nm (1-m). 
The joining of zero lines in (c) results in the formation of a new lasing mode (new zero line crossing 
is encircled in white). The inset in (c) is an enlargement of the mode 17 and new mode solutions. 
In (d), the phase singularity at the new mode has a topological charge of +1, opposite to that 
of mode 17. The real and imaginary zero lines pull apart in (f) so that the mode 17 and new 
mode solutions are nearly identical. The phase map in (g) reveals the existence of the two phase 
singularities. The lines completely separate in (h) resulting in the disappearance of mode 17 and 
the new mode. The phase map in (i) confirms the disappearance of the two modes. This process 
reverses itself in (j-m) yielding behavioral symmetry around Iq = 14472 nm. 
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FIG. 9: (Color online) (a) Standing/traveling wave ratio Ast( x ) of the new lasing mode (red) 
and lasing mode 17 (black) for Iq = 14553 nm (solid lines) and Iq = 14523 nm (dotted lines). The 
potential profile |VT"(x)| 2 (black dashed line) of the dielectric function for this wavelength is overlaid 
in both (a) and (b) and major potential barriers are marked © through ©. The ratio A$t{x) of 
the new lasing mode and lasing mode 17 become more similar an converge on each other as Iq 
reduces. Reducing Iq further causes these two lasing modes to first disappear then reappear as the 
process reverses itself, (b) Standing/traveling wave ratio A$t{x) of the new lasing mode (red) and 
lasing mode 17 (black) after they have reappeared for Iq = 14284 nm. The ratios are similar to 
the ratios for Iq = 14553 nm in (a). The vertical black solid line marks the gain edge. The ratios 
of the modes diverge now as Iq is reduced. 



singularity in Fig. |8[d). The new mode is close to mode 17 in the (k, n») plane and its 
phase singularity has the opposite topological charge as seen in Fig. |8](d). As l G decreases 
further, the joined zero lines forming mode 17 and the new mode pull apart. This causes the 
two solutions to approach each other in the (k, n*) plane, i.e., the frequency and threshold 
of mode 17 increase while the frequency and threshold of the new mode decrease. In Figs. 
8^f) and (g), the solutions are so close that they are nearly identical, yet they still represent 
two separate solutions. Further decreasing Iq makes the solutions identical. The zero lines 
then separate and the phase singularities of opposite charge annihilate each other in Figs. 
8pi) and (i). This results in the disappearance of mode 17 and the new mode. The process 
then reverses itself as Iq is decreased further [Figs. [8]j)-(m)] yielding the reappearance of 
mode 17 and the new mode and their subsequent separation in the (k, rij) plane. This is 
the aforementioned behavioral symmetry around Iq = 14472 nm. 
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Examining the standing/traveling wave ratio of lasing mode 17 and the new lasing mode 
together with the potential profile |lV(x)| 2 offers some insight of mode annihilation and 
reappearance in real space. Figure [9] shows the ratio Ast{%) for the new mode and mode 
17 along with ^(x)) 2 . The potential profile is very similar for the new mode and mode 
17 since their wavelengths are very close. There are four major potential barriers at the 
mode 17 wavelength (A = 598 nm) for x < 15000 nm. This is the spatial region associated 
with the gain distributions in Fig. [8] where Iq is always smaller than 15000 nm. Figure 
[9] shows them at: x = © 927 nm, © 5200 nm, © 8700 nm, and © 14519 nm. Due to 
oscillations, the centers of barriers © and © are less well defined. The right edge of the gain 
region at Ig = 14553 nm is located just to the right of barrier ©. For Iq = 14523 nm, the 
right edge of the gain region nears the maximum of barrier ©. Figure |9](a) shows that for 
Ig = 14553 nm, the SWC of the new mode is between barrier © and barrier ©. The SWC 
of mode 17 is in the middle of the gain region at x = 5300 nm and its SWC is between 
barrier © and barrier ®. Before disappearing, the modes approach each other in the (k, rij) 
plane, eventually merge, and their intensity distributions become identical (as evidenced by 
the trend of their standing/traveling wave ratios). As Ig is further reduced and the modes 
reappear, the behavior of the modes' ratios Ast{%) (or equivalently, intensity distributions) 
reverses itself as expected from the behavioral symmetry shown in Fig. [8j At Ig = 14284 
nm, the right edge of the gain region has passed barrier © and Fig. |9|b) [with a different 
horizontal scale than Fig. [9]^a)] shows the SWC of the new mode is in roughly the same 
location as it was for Iq = 14553 nm. The SWC of mode 17 is also in roughly the same 
location as it was for Ig = 14553 nm. 

The appearance of new lasing modes is unanticipated. In the passive system, the number 
of standing wave peaks for quasi modes increases incrementally by 1, e.g., quasi mode 17 
has 82 peaks and quasi mode 18 has 83 peaks. Lasing modes 17 and 18 behave the same 
way. How exactly does a new lasing mode fit into this scheme? Though closer in frequency 
and threshold to lasing mode 17, counting the total number of standing wave peaks of the 
new lasing mode yields the same number as for lasing mode 18. However, the new lasing 
mode is somewhat compressed in the gain region having one more peak than lasing mode 
18. It is decompressed in the region without gain having one less peak than lasing mode 18. 

Comparing the decompositions of the lasing modes in terms of quasi modes helps reveal 



the character of the new lasing mode. Figure 10 shows the decomposition of the new lasing 
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FIG. 10: (Color online) Decomposition at Iq = 14284 nm of lasing mode 17 (black circles), lasing 
mode 18 (blue open diamonds), and the new lasing mode (red crosses) in terms of the quasi modes 
of the passive system. Lasing modes 17 and 18 are mostly composed of their respective quasi modes 
while the new mode is dominated by a mixture of both quasi mode 17 and 18. The inset shows the 
decomposition coefficients of outlying quasi modes for lasing mode 17 (black line), lasing mode 18 
(blue line), and the new lasing mode (red line). 

mode together with the decomposition of lasing modes 17 and 18 at Iq = 14284 nm. The 
new mode has a slightly larger coefficient amplitude associated with quasi mode 17 than 
quasi mode 18, but the two amplitudes are nearly equal. We found that as mode 17 and the 
new mode solutions approach each other by varying Iq, their coefficient distributions also 
approach each other until becoming equal as expected from Figs. [8] and [9] 

IV. CONCLUSION 

We have demonstrated the characteristics of lasing modes to be strongly influenced by 
nonuniformity in the spatial gain distribution in ID random structures. While the entire 
structure plays the dominant role in determining the frequency of the lasing modes, the gain 
distribution mostly determines the lasing thresholds and spatial distributions of intensity. 
The gain distribution also appears to be solely responsible for the creation of new lasing 
modes. We have verified the existence of new lasing modes in numerous random structures 
as well as dielectric slabs of uniform refractive index. A more thorough investigation of the 
latter will be described in a future work. All of these changes caused by nonuniform gain 
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take place without the influence of nonlinear interaction between the field and gain medium. 
Our conclusion is that nonuniformity of the gain distribution alone is responsible for the 
complicated behavior observed here. 

By decomposing the lasing modes in terms of a set of quasi modes of the passive system, 
we illustrated how the lasing modes change. The contribution of a quasi mode to a lasing 
mode was seen to depend mostly on its proximity in frequency k and the spatial distribution 
of gain. The more the gain changed from uniformity, the greater the mixing in of neighboring 
quasi modes. Thus, great care must be taken even close to the lasing threshold when using 
the properties of quasi modes to predict characteristics of lasing modes in weakly scattering 
systems with nonuniform gain or local pumping. 

The change of intensity distributions of lasing modes as the size of the gain region is 
varied appears to be general. With reduction of the size of the gain region, the peak of 
the standing/traveling wave ratio A ST (x), or the standing wave center (SWC) of the mode, 
moves to stay within the gain region. Modes with low thresholds have a SWC near the middle 
of the gain region while high threshold modes have a SWC near the edge of the gain region. 
Changing the gain distribution thus changes the intensity distributions of lasing modes. The 
exact modal distributions, however, appear correlated with the potential profile. In the cases 
studied here, the new lasing mode and lasing mode 17 lay in between two large potential 
barriers. Decreasing the size of the gain region brought the intensity distributions closer 
together until they disappeared. These changes took place by varying the edge of the gain 
region only hundreds of nanometers. Thus, even a slight change in the gain distribution 
may have drastic consequences for lasing modes. 
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APPENDIX A: LINEAR GAIN MODEL 



In this appendix, we describe the model used to simulate linear gain in a ID system. The 
gain is linear in the sense that it does not depend on the electromagnetic field intensity. The 
lasing solutions ty(x) must satisfy the time-independent wave equation 

d 2 



— + e{x ) u)k 2 



= 0, (Al) 
with a complex frequency-dependent dielectric function 

e(x,u) = e r (x) + x 9 {x,u), (A2) 

where e r (x) = n 2 (x) is the dielectric function of the non-resonant background material. The 
frequency dependence of e r (x) is negligible. cj), corresponding to the susceptibility of 

the resonant material, is given by 

A e N A (x) 

Xg\X,U ) = — 2 2 ~A ' A3 

where A e is a material-dependent constant, Na(x) is the spatially dependent density of 
atoms, u a is the atomic transition frequency, and Au a is the spectral linewidth of the 



atomic resonance. Equation (A3) may be simplified by assuming the frequencies of interest 



uj are within a few linewidths of the atomic frequency u a , i.e., u 2 — u 2 = {u + uo a ) (uj — uj a ) 



2u) a (u) — u a ). Equation (A3) then reduces to 



uj a Au a [l + 2i(u - u a )/Au} a \ 
The frequency-dependent index of refraction is 

n(x,uj) =^/e(x,u) = Jejx) + Xg(x,u) 

—n r (x,uj) + irii(x,ui), (A5) 

which may then be implemented in the transfer matrix method. At this point, let us note 
that only 2 steps are needed to convert this classical electron oscillator model to real atomic 



transitions |29j. First, the radiative decay rate 711 may be substituted in to Eq. (A4) in 
place of a few constants. Second, and more importantly, real quantum transitions induce 
a response proportional to the population difference density AN a- Thus, Na{x) should be 
replaced by AN a, the difference in population between the lower and upper energy levels. 
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Linear gain independent of u is obtained by working in the limit uj — uo a <C Acj a , yielding 



.A e &N A (x) 

l 7 



(A6) 



a purely imaginary susceptibility. We can make the definition Xg( x ) = i e i( x ): where ti{x) is 
the imaginary part of e(x). Note that e(x) may include absorption [e* > 0] or gain [e$ < 0]. 
We shall only consider gain here. The complex frequency-independent dielectric function 
now yields a frequency-independent index of refraction n(x) = n r (x) +irii(x) which may be 
expressed explicitly as 



n r (x) 



Ui(x) 



n[x) 



V2 



—n(x) 



1 + 



e?(x) 



n 4 (x) 



1/2 



V2 



n 4 (x) 



1 



1/2 



(A7) 



Furthermore, in the main text, we assume ni to be spatially independent. Thus, by solving 
for n r (x) in terms of n(x) and n^, the index of refraction used throughout this paper is given 
by 



n(x) —n r (x) + ini 

= \Jn 2 (x) + nf + ini. 



(A8) 



APPENDIX B: STANDING WAVE AND TRAVELING WAVE COMPONENTS 
OF THE TOTAL FIELD 

In this appendix, we describe the method that enables one to define a standing wave 
component and a traveling wave component of the field at each point a; of a ID system. 
For an open structure without gain, the field reads 

if)(x) = p(x) exp[in(x)kx] + q(x) exp[—in(x)kx], (Bl) 

where k is the complex wavevector and n(x) is the index of refraction, the value of which 
alternates between n(x) — n\ > 1 in dielectric layers and n(x) = n 2 = 1 in air gaps. For 
structures with gain, the field reads 

ty(x) = p(x) exp[ih(x)kx] + q(x) exp[—ih(x)kx], (B2) 
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where n(x) = n(x) + m< is the complex index of refraction. We rewrite both equations in 
the single form 

E(x) = p(x) exp[iK(x)x] + q(x) exp[—iK(x)x], (B3) 

where K(x) = K r (x) + iK{(x) and E(x) may be either ip(x) or ty(x). 

For now, we will consider the field within a single layer in order to simplify the notation. 
The following results will be valid within any layer. Since within a layer, the coefficients 



p(x), q(x) and the wavevector K(x) do not depend on x, we rewrite Eq. (B3) as 



E(x) = pexp[iKx] + qexp[— iKx]. (B4) 

The complex amplitudes p and q of the right-going and left-going fields, respectively, can be 
written as p = P exp[i(p] and q = Qexp[i(p} where P and Q are the real amplitudes which 
can be chosen positive. The field becomes 

E(x) =Pexp[— Kix] exp[i(K r x + cp)] 

+ Q exp[Kix] exp[— i(K r x — 0)] 
=II(x) exp[i(K r x + ip)} 

+ Q{x) exp[-i(K r x -(f))}, (B5) 

where U(x) = Pexp[— K^x] and Q(x) = Q exp[Kix]. Introducing the global phase $ = 
[ip + 0]/2 and the difference A = [ip — <p}/2, the field reads 

E(x) =exp[2$]{Il(x) exp[i(K r x + A)] 

+ Q(x) exp[—i(K T x + A)]}. (B6) 

Within a single layer, we can set $ = so that the field becomes 

E(x) =n(x) exp[i(K r x + A)] + 0(x) exp[-i(K r x + A)] 

=E iR \x) + E^(x), (B7) 

where E^(x) and E^ L \x) are the right-going and left-going waves, respectively. 
We can build a standing wave component with E^ R \x) as 

E^ s \x) =EW(x) + [E( R \x)]* 

=2U{x) cos[K r x + A] (B8) 
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and define the traveling wave component as the remaining part of the total field 



E {T \x) =E(x)-E {s \x) 

=E iL \x) - \E( R \x)}* 

= [6(x) - U(x)] exp[-i(K r x + A)]. 



(B9) 



Hence, 2U(x) and — are the amplitudes of the standing wave and traveling 

wave components, respectively. It is also possible to build a standing wave component with 
E^ L \x) as 



E^ s \x) =E i ~ L \x) + [E i - L \x)}* 
=290) cos[K r x + A] 



(BIO) 



so that the traveling wave component reads 



E (T \x) =E(x) - E (s \x) 

( x ) - [E^ L \x)}* 
= [n(x) - Q(x)] exp[i{K r x + A)]. 



(Bll) 



Comparing both ways of resolving the total field into its two components, we see that 



in Eq. (B9) the traveling wave component is a left-going wave while in Eq. (Bll) it is a 



right-going wave. Hence, if in the expression of the field in Eq. (B7), the prevailing wave is 



the right-going wave H(x) exp[i(K r x + <p)\ (i.e., H(x) > 0(x)), we choose the standing and 



traveling wave components of Eqs. (BIO) and (Bll). In the opposite case of H(x) < 0(x), 



we choose the standing and traveling wave components of Eqs. (B8) and (B9). 



Let us note that the imaginary part of the total field E(x) is given in both cases by 



lm[E(x)] = [U(x) - Q(x)] sm[K r x + A]. 



(B12) 



As expected, the presence of a traveling wave component, i.e., |n(x) — 0(x)| 7^ 0, makes 
E(x) become complex instead of being real for a pure standing wave. 
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